function z=f31( x, y )

	z = ( 3 * ( x ^ 2 ) + 2 * x ) / ( 8* y );

endfunction
